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Oh 1 Abstract 

q 

^ We extend the theory of Talbot revivals for planar or rectangular geometry to the 

• i-H . case of cylindrical waveguides. We derive a list of conditions that are necessary 

to obtain revivals in cylindrical waveguides. A phase space approach based on the 

_C Wigner and the Kirkwood-Rihaczek functions provides a pictorial representation of 



TM modes interference associated with the Talbot effect. 
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1 Introduction 



O 

Although the Talbot effect was discovered firstly in the beginning of nineteen 
century (1836) and then rediscovered many times in different systems, it has 
never been a widely known phenomenon. In quantum systems it is often re- 
ferred as "self-imaging" or "quantum revivals". Even in the optical domain 
not always Talbot name is mentioned when this phenomenon is described 
and discussed. However, it was W. H. F. Talbot, who firstly observed that 
monochromatic light passing through periodic grating at a certain distance 
from the grating forms its ideal image and consecutively at integer multiples 
of this distance similar images are reproduced. He also demonstrated that 
when white light is used, images of the grating of different colors are formed 
at different distances and that is why we believe that all "self-imaging" effects 
should be referred as Talbot effects. 

After the original paper by Talbot [1] , followed by the work of Lord Rayleigh 
[2], and then by a series of papers of Wolfke [3], many papers have been writ- 
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ten about this subject. A comprehensive description of the Talbot effect and 
its rediscoveries in classical optics can be found in Patorski [4]. Similar his- 
torical review and a detailed description of "quantum revivals" can be found 
in [5]. In the last years of the XX century fractional aspects of the Talbot 
effect attracted considerable interest [6,7,8,9]. These effects were studied both 
in optical and quantum mechanical domains. As regards Talbot effect in op- 
tical waveguides most significant are works of Ulrich [11], who investigated 
waveguides of planar and ribbon geometries. 

The literature concerning the Talbot effect in planar or rectangular waveguide 
geometries is quite reach (even some US patents for applications of the Talbot 
effect in those systems exist), but the much more practical cylindrical geom- 
etry was not taken into account in optical studies. We present in this paper a 
comprehensive theoretical study of the Talbot effect in cylindrical waveguides. 
This should fill the gap in the available descriptions of self-imagining phenom- 
ena in various waveguides. There is a limited number of papers mentioning 
cylindrical geometry in quantum revivals [12], but the results presented in this 
paper go beyond the one obtained so far. 

The paper is organized as follows. In Section 2 we present general assumptions 
from which our study starts. We use a formal similarity between field prop- 
agation in waveguides and the dynamics of wave packets in potential wells 
in a given geometry. Using this analogy phase space Wigner and Kirkwood- 
Rihaczek functions are used. The phase space functions provide a pictorial 
description of interference effects, that are the basis of the Talbot revivals. In 
Section 3, the case of dielectric fibers and the possibility of revivals in this most 
practical for possible application system is studied. In Section 4, solutions of 
the wave equation in cylindrical mirror waveguides are analyzed and require- 
ments that have to be fulfilled to obtain revivals are derived and imposed. 
A detailed study of approximations used is included. Finally, a summary of 
the results is given with a short paragraph devoted to the presentation of the 
applied method. 



2 Talbot effect in phase space 

2.1 General assumptions 



The reason why the Talbot effect appears both in optical and quantum me- 
chanical systems, mathematically can be summarized very briefly: The Helmholtz 
equation is common for classical electrodynamics and quantum mechanics. 
The dependence of the electromagnetic field in the direction of field propaga- 
tion can be regarded as an analogue of time dependence of a wave packet in 



quantum mechanics. In the regime where paraxial approximation is justified 
this analogy is especially clear. 

This paper is focused on Talbot revivals of an initial field in cylindrical waveg- 
uides which is purely an optical phenomenon but, nevertheless, as we shall see 
in the next paragraphs, its most simple explanation can be given referring to 
quantum mechanical concept of phase space distributions. The key question 
we shall pose and answer is whether the Talbot effect in cylindrical waveguides 
exists and can be used in practice. 

We assume that harmonic, monochromatic plane waves propagate through the 
waveguide which symmetry axis was chosen as the z direction of the system. 
Inserting the fields 

E(x, y, z, t) = E(x, y)e ±lkz ^\ B(x, y, z, t) = B(x, y)e ±ikz -^ , (1) 

into the Maxwell equations, we obtain the two-dimensional Helmholtz equa- 
tion: 
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= , where Y = ^e— - k 2 , Vl = V 2 - 



c 2 ~ dz 2 



The propagating constants /c, corresponding to specific modes are to be derived 
from appropriate boundary conditions (see, e.g. [13]). 



2.2 Phase space distributions 



As we have already mentioned the basis of Talbot effect, i.e. interference, is 
especially clearly seen in the phase space description. The quasi-distribution 
functions used for the study of electromagnetic fields, especially pulses, usu- 
ally work in the time-frequency domain. Here, a totally different approach was 
chosen as all the fields taken into consideration are assumed to be monochro- 
matic and of a harmonic time-dependence. Thus, we shall use a typical phase 
space known from mechanics constructed from the position and momentum 
variables. 

Some complications result from the fact that the electromagnetic field requires 
a vector description. However, in many cases (e.g. mirror waveguides) solutions 
of the Maxwell equations are divided into TM and TE modes, i.e. are of 
type entirely determined by the E z or B z field component, respectively. Thus, 
we can treat E z (B z ) component, that fully describes the field, as a single 



scalar function corresponding to the quantum mechanical wave function of a 
potential well problem with analogous geometry. 

To present problems of revivals is phase space we shall use the following two 
quasi-distribution functions: the Wigner function [14,15], 



W*(x,p) 



(2vr) c 



#*(£ + f/2) e iR V(x - |*/2) d d £ 



(2) 



which is the most commonly known phase space distribution, and the Kirkwood- 
Rihaczek (K-R) function [16,17]: 



K(x,p) 
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that is very convenient for calculations. Subscript d in the definitions above de- 
notes number of spatial dimensions of the system, whereas ^ is a function that 
characterizes a system (e.g. wave function in quantum mechanical applications, 
or E z field component for TM modes in waveguides, etc.) Because for 3 dimen- 
sional systems phase space distributions are 6 dimensional, graphically only se- 
lected cross-sections can be presented. This arbitrary choice is made somehow 
easier for solutions of the form of Eq. (1), when the z-dependence of the field 
separates from the transversal components. Figure 1 shows examples of the 
"transversal" Wigner and K-R functions calculated for E z = Ja(jir) — Jo(J2r) , 
which corresponds to superposition of TM i and TM 2 modes in cylindrical 
mirror waveguide: j\, j'2 denote the first and the second zero of Bessel function 
Jo(u) and r = J x ~ty , where a is a waveguide radius. 



a) 





Fig. 1. The y = and p y = cross-sections of "transversal" quasi-distributions 
corresponding to E z = Jo(jir) — Jo(j2?"): a ) cross-section of the Wigner function; 
b) cross-section of the real part of the K-R function. 



As we have already mentioned in the case of solutions of Maxwell equations 
having the plane wave form we consider here, the ^-dependence of the field sep- 
arates from the transversal components. Thus, when describing Talbot effect, 



we can simply "forget" the transversal dependence of the field and concentrate 
on the z-p z cross-section of the phase space distribution which is essential for 
the Talbot revivals. Integrals corresponding to the exp(— ikiz) factors are quite 
elementary and the cross-sections of quasi-distribution functions for given x, 
y, p x , p y , i.e. set the transversal position and momentum components, are easy 
to obtain. 

The Wigner function for a superposition of two plane waves exp(— ik\z) + 
exp(— ik 2 z) is given by 



W 2 (z, Pz ) ~ 5( Pz - h) + 5( Pz -k 2 ) + 25 (p z - h±tl\ cos ( z (k 2 - h)). (4) 
The real part of the K-R distribution for such a superposition takes the form 



Re[K2(z,p z )]~5(p z -k 1 )+5(p z -k2) + [5(p z -ki)+5(p z -k2)]cos[z(k 2 -k 1 )]. (5) 
When there are N superposed waves, the following sums are obtained: 



i=l l=i<j=N ^ ' 



W N (z,p z )~y ]S(p z -ki) + 2 Jj H Pz - ^r^ ) cos[z(k t - kj)} (6) 

i<j= 



for the Wigner function and, for the real part of the K-R function, 



Re[KN(z,pz)]~J26( Pz -k j )(l+ V coa[z(ki-kj)]), (7) 
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with appropriate coefficients. In both cases the whole z dependence is inserted 
into interference cos[z(ki — kj)] terms. Initially, at z = 0, all this cosines are equal 
to 1. The further z dependence is guided by (ki — kj) factors. When all these 
(hi — kj) factors are commensurable with each other, perfect regular revivals 
are obtained at such z rev for which all cos[z rev (ki — kj)] are simultaneously equal 
to 1 again. 

Obviously, all superpositions of just two different modes, like Eqs. (4), (5), will 
revive perfectly at multiples of , 2 * , no matter what kind of waveguide we 
consider. But, when there is more superposed modes, commensurability of all 
possible (ki — kj) factors is needed to obtain perfect revivals. As the propagat- 
ing constants ki vividly depend on the type of waveguide and its parameters, 
such a commensurability is rather an exceptional then a typical case. Even 
dealing with highly symmetric problems like mirror waveguides of planar or 



square cross- sections we have to keep in mind that ki = Jk$ — jf which means 
that commensurability of 7j's (in these waveguides 7j's are commensurable) 



is not automatically followed by commensurability of propagating constants 
hi. Only when linear approximation of Jk\ — jf holds, commensurability of 
7i's is sufficient and that is why we shall often limit ourselves to the lowest 
from propagating modes. A definite advantage of working within the range 
where the linear approximation of square root holds for systems having 7j's 
proportional to each other is the fact that for a given wavelength the Tal- 
bot distance is settled, it does not depend on superposed modes. Otherwise, 
different modes superpositions shall revive at different distances. 



All the features characteristic for the Talbot effect can be clearly explained by 
looking at interference cos[z(ki — kj)} terms, their getting in and out of phase. 
The phase space description brings us in a natural way to this simple idea and 
indicates how fundamental this concept is. 



3 Talbot effect in dielectric waveguides 



Cylindrical dielectric waveguides are called optical fibers. We shall analyze 
only the step-index fibers, i.e. the fibers with constant refractive indexes in 
the core and the cladding, which are entirely characterized by radii of core and 
cladding and their reflective indexes n\ = y/JIiei and n 2 = v //i2^2- Assuming 
that the cladding and the core differ only by dielectric constants (magnetic 
permeabilities fii = 1x2) the standard boundary conditions lead to the following 
equation [18,19] 
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where a denotes the core radius, r = cry , k — af3 , 7 = ^iEi^t — k an d 
1 = k 2 — fi>2£2^r- Although Eq. (8) has a quite nice regular form there is no 
way to solve it analytically. In the simplest case when the field has no azimuthal 
dependence, i.e. v — 0, its solutions can be divided into TE and TM type, as 
for the mirror waveguides. Then a graphical picture gives a clear representation 
of solutions similarly to the case of finite potential wells in quantum mechanics. 
However, we have to keep in mind that the (^—independence, corresponding to 
v — 0, is not a typical case. General solutions of Eq. (8) are if dependent and, 
actually, the lowest propagating mode in step-index dielectric fiber is obtained 
for v — 1. Thus, let us start our analysis from general solutions. 



3.1 General solutions in step index fibers 



The lowest propagating mode in a cylindrical dielectric waveguide is always an 
HEu mode (in this notation HE means that field H z dominates over E z field, 
while for EH modes the E z field dominates). Single mode fibers are of a great 
practical importance, but in this study we are not interested in them because a 
single mode propagates without a change in its transverse distribution. Next 
modes are TEqi, TMqi, and HE21. They appear almost simultaneously as 
their propagating constants are nearly the same. 

Let us firstly consider a fiber in which there are only these four propagating 
modes (e.g. A = 1550nm, n\ = 1.46, n 2 = 1.45, a = 4.5/xm). We shall refer to 
this situation as to the "limit of small number of modes" . Numerical solutions 
of Eq. (8) obtained for these parameters are: a'jHEu = 1.79268, o^te 01 = 
2.75973, Gr/TMoi — 2.76234, a^HE 2 i — 2.76342. While analyzing problem of 
revivals of superpositions of HEu, TE i, TM 01 , and HE 2 ± modes, one has 
to realize that although four modes are superposed, it is a special case when 
three propagating constants (corresponding to TE 01 , TM m , and HE 21 modes) 
are really close to each other. As we have already mentioned, superpositions 
of arbitrary two modes with propagating constants fci, k 2 shall revive at all 
integer multiples of zt = k 37 ■ It is easy to calculate that initial images 
constructed from superposition of pairs (HEu, TyE 01 ), or (HEu, TM 01 ) or 
(HEn, HE 2 i) shall revive at multiples of z\ = 3,40787mm, z\ = 3.3967mm, 
z\ = 3.39214mm, respectively. Those values are so close to each other that 
images constructed from all four modes shall revive at a mean value which 
is zt = 3.39891mm. In numerical simulations of field propagation one can 
observe the first, second and even 20th or 50th Talbot revival. Obviously, 
higher revivals are less accurate, and after several tenths of faithful revivals 
they get out of phase - but then, after some propagating distance they get in 
phase again, and again some faithful revivals are obtained. Then the situations 
of getting in and out of phase appropriate cosines repeats almost cyclicly. 

Plots of initial light intensity corresponding to superposition of TM i—TE i — 
HEu — 2HE 2 i modes, its first Talbot revival and a characteristic inverted 
image at the half of the Talbot distance are presented in Figure 2. Infidelities 
of subsequent revivals, can be calculated using standard measure 



||J(0)-J(z)|| ( Cd<pfidpp\I(p,<p;Q)-I(p,<p;z)\* y 

MZ) ll'(o)ll lv f?*pSJdp P \i(p,<r,o)\* ) ' w 

to which we refer as to the infidelity as its value increases with increasing 
deviation of I(z) from original 1(0) and only when / = the copy is perfect. 
For superposition of HEu — TE i — TM m — 1HE 2 \ modes the infidelities 
of the successive revivals are quite low (/ = 0.0000276, /o = 0.0029825, 



/o = 0.0260536, f = 0.0945819 correspond to the 1st, 10th, 30th, and 60th 
revival, respectively) and in such a four-mode fiber infidelities of revivals for 
every initial field distribution would be of this order. 



/o =0.0000276 




Fig. 2. Superposition of TMqi — TE$\ — HE\\ — 2HE21 modes: a) the original 
intensity; b) intensity at the half of the Talbot distance; c) 1st Talbot revival. 

It is seen that in the limit of small number of modes Talbot revivals in dielectric 
fibers can be obtained. Change of fibers parameters will result in the change of 
Talbot distance, but it is quite obvious that for a small number of propagating 
modes Talbot revivals can be obtained. Just as we have seen using phase space 
representation - appropriate number of cosines have to get in phase to obtain 
the revivals. 

This simple and quite intuitive method of looking for the Talbot distance starts 
to be more complicated when the number of propagating modes increases. 
Then, different methods are have to be applied to calculate revival distance 
as we shall see in the next section on the example of TE 0i modes. 



3.2 The ip -independent case (v = 0) 



The ^-independent solutions of Eq. (8) can be divided into TE and TM 
modes. For TE modes the following relation is obtained: 



Ux(r) i k^k) 
rjo(r) kK (k) 



0. 



(10) 



Because the parameters Y and k are (by definition) correlated, 

n 2 , i 2 

r + k — — — (717 — ni) := V , 



(11) 



we obtain a set of equations that can be solved graphically or numerically. 
It is convenient to introduce a normalized frequency parameter V, Eq. (11), 
and depict for example — ^ j44 and ^ A - 1 [^| as a function of Y for appropriate 
values of V. 



For TM modes instead of Eq. (10) we would have 

Ui(r) eil gipc) 

rJo(r) s 2 kK (k) 



0. (12) 



Graphical solutions of Eq. (12) are similar to those of Eq. (10), the only 
difference is that the Macdonald part of the plot kI ^ ^\ is modified by a fixed 
factor — , which is usually close to 1. In practice this means that TE 0n and 
TM Qn modes have nearly the same propagation constants and they will tend 
to appear simultaneously. 

In the previous subsection we were dealing with the small number of modes, 
so now let us focus on a limit of large number of propagating modes (large 
frequency parameter V). In the limit V — ► oo, solutions Gry n would be given 
by zeros of the Bessel J\(p) function, and propagating constants would cor- 
respond exactly to those obtained for TE$ n modes in the mirror waveguides. 
This case would be widely discussed in the next Section, thus, now we shall 
take into consideration only the finite values of V. Analyzing the graphical 
representation of Eq. (10) for different values of frequency parameter V one 
finds that with the increase V the Macdonald part of the plot, kI ^ ^K , starts 
to be parallel to T axis for the increasing range of T n 's. Moreover, it is also 
getting closer to this axis as for T — > value of K ^\ — y- For large V 
and low mode numbers solutions would be of very regular form: the first one 
corresponding to TE 0i mode will be given by, say T , and the approximate 
formula for the next solutions would be T n = Tq + amr. Obviously the larger 
V, the more accurate this formula is, and it works well only for modes having 
numbers low in comparison to the total number of modes. Thus, in order to 
obtain revivals we shall use only a few percent of the lowest from propagating 
modes for constructing initial images, namely those modes for which the linear 
approximation of the square root is sufficient. Sometimes, however, it is easier 
to omit such analytical approximations and to simply calculate numerically 
infidelities of intensity distribution as a function of z and look for the minima 
of this function. 



3.2.1 Examples of revivals 

In this Section we consider only the (^—independent fields that are fully charac- 
terized by their cross-section along the radius. Thus, the figures presented shall 
show cross-sections of the light intensity versus the normalized distance from 
a fiber center p/a. All the examples were calculated numerically for a quite 
thick fiber (a = lmm) with refractive indexes of the core and cladding equal to 
rii = 1.47, i%2 = 1.45, respectively and a wavelength of A = 850nm. For these 
values of A, rii, and n 2 the frequency parameter equals V= 1786.35, which 
correspond to more then 550 propagating TE modes, and 7 = 3.82956/a. 



a) symmetric superposition of TEqi + TE02 + TE03 + TEqa modes 
Figure 3 presents the light intensity corresponding to the superposition TE m + 
TE 02 + TE 03 + TE 04: at z = 0, z = z t = 27.7036m, z = 10z t , z = 30z t , z = 60z t , 
and z = 100z t . The original intensity is depicted in black, the revivals in blue. 
Above the plots of revivals their infidelities are depicted. The value of the 
Talbot distance z t was determined numerically by finding minimum of the 
infidelity function, Eq. (9). It is seen that although revivals are not perfect they 
are certainly faithful enough even at distances of lkm. For comparison, Figure 
4 presents examples of the light intensities at distances between revivals: it is 
clear that the intensities at multiples of the Talbot distance differ significantly 
from typical intensity distribution during propagation. 



0.00085 
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Fig. 3. Light intensity for the symmetric superposition of TEq u for n G {1, ...,4}. 
The original intensity distribution are plotted in black, the intensity at the Talbot 
distance zt = 27.7038m and its multiples in blue, a) the original intensity; b) 1st 
Talbot revival; c) 10th Talbot revival; d) 30th Talbot revival; e) 60th Talbot revival; 
f) 100th Talbot revival. 
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Fig. 4. The same superposition as in Fig. 3. The original intensity distribution is 
plotted in black, in orange are plotted intensities at a) z = lm, b) z = 5m, c) 
z = 10m. 



b) Gaussian "25" 

Figure 5 presents the revivals of the initial intensity having a Gaussian radial 
distribution. The Gaussian function was obtained from a superposition of 
the first 25 TE Q _ modes. The numerically calculated Talbot distance is z t = 



10 



27.7035m and the infidelities of revivals (depicted above every plot) are quite 
low. It is clear that revivals of intensities of a given shape can be observed. 
Obviously, this result is true only on the assumption that other modes (e.g. 
modes depending on ip) do not contribute to the initial image. However, the 



4*10" 



/o =0.00082 




0.8 
0.6 

0.4 
0.2 



0.4 0.6 0.8 1 b) 
/o =0.00730 



0.2 0.4 0.6 0.8 1 
/o =0.02798 



0.2 0.4 0.6 0.8 1 
/o =0.07252 




d) ■ • • ■ e) 

Fig. 5. Light intensity for the superposition of TEq h , n £ {1, ...,25}, providing a 
Gaussian distribution. The original intensity distribution are plotted in black, the 
intensity at the Talbot distance zt = 27.7035 and its multiples in blue, a) original 
intensity; b) 1st Talbot revival; c) 10th Talbot revival; d) 30th Talbot revival; e) 
60th Talbot revival; f) 100th Talbot revival. 

fidelities of revivals depend significantly on the effective number of terms of 
the Bessel- Fourier (BF) series contributing to the initial image. Although in 
the example presented above we have taken first 25 terms of the BF series, 
only first 9 coefficients were larger then 1/1000 and only first 4 were larger 
then 1/100. The coefficients from 11th to 25 were of the order of 5/10000. 
If we would like to propagate, say, a more slim Gaussian these proportions 
would be different. Is is worth noting that in the case of dielectric waveguides 
this effective number of contributing modes effect not only the fidelities of 
revivals, but also the optimal Talbot distance, which we will clearly see in the 
comparison with the next example. 



c) symmetric superposition of TE m + TE 02 + • • • + TE 02 $ modes 
This example of revivals in a thick fiber illustrates how the effective number 
of contributing modes might modify a Talbot distance. As the initial intensity 
we take the one corresponding to the symmetric superposition of TEqi + 
TE 02 + - ■ - + TEq25 modes. Figure 6 presents the initial intensity, its 1st Talbot 
revival and then 10th, 30th, 60th, and 100th Talbot revival. The numerically 
calculated "optimal" Talbot distance is equal in this case z t = 27.7013m, 
which is slightly smaller then in the case presented in examples a), b). 

We have already stressed the fact that different initial images have different 
Talbot distances can be explained by the effective number of modes that are 
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Fig. 6. Light intensity for the symmetric superposition of TE§ n modes for 
n € {1, ..., 25}. The original intensity distribution are plotted in black, the intensity 
at the Talbot distance zt = 27.7013m and its multiples in blue, a) the original in- 
tensity; b) 1st Talbot revival; c) 10th Talbot revival; d) 30th Talbot revival; e) 60th 
Talbot revival; f) 100th Talbot revival. 



superposed. It is, again, a consequence of importance of the quality of approx- 
imations used. In our simplified analysis we have assumed that 7 n ~ 70 + mc. 



On the one hand this is true only for n small enough for 



KjJK,) 
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parallel to the 07 axis, on the other distance between neighboring zeros of 
Bessel functions is ~ ir only for higher modes. Thus, it is obvious that if only 
first few modes contribute to the initial image (or their contribution domi- 
nates) the real Talbot distance would be slightly different then if there were 
more contributing modes. Differences in optimal Talbot distance z t for su- 
perpositions presented in examples a), b), c) can be evaluated explicitly by 
calculating the infidelities of revivals. They start to be important at 20th, or 
50th revival when the initial 2mm difference in optimal Talbot distances re- 
sults in 5 or 10cm divergence from the distance of optimal revival. As long as 
we are interested only at the first Talbot revival we can take an average z t and 
the infidelities of revivals of different initial images at z t should not be larger 
then 1/100 which is sufficient for most applications. Obviously, when we are 
interested in revivals at larger distances all the initial images can be divided 
in classes having the same "effective number of contributing modes" and the 
revivals at multiples of the corresponding "optimal Talbot distances" would 
be obtained for all the images within the class. 



We have shown numerical simulations indicating that the Talbot revivals of 
initial images constructed from TE 0n modes can be obtained in optical fibers. 
Although examples a) - c) present superposition of TE modes, one can note 
that for TM 0n modes revivals should be even more faithful, as — < 1 and thus 
7„ are even closer to zeros of J\{p) function then it is for TE modes, which 
brings us to the case of mirror waveguides. 
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4 Talbot effect in mirror waveguides 



Although we present it in the last Section, a model of ideal mirror waveguides 
is very useful for preliminary calculations. It is analytically soluble for systems 
of standard geometries because boundary conditions are quite simple: normal 
component of B and tangential component of E have to vanish at the boundary 
mirror surface. 



4-1 Planar mirror waveguides 



In the elementary case of planar mirror waveguide the propagation constant 
for nth mode is given by k\ = k$ — 1L jh, where d denotes a separation distance 
between mirrors plates [20,18]. It is clear that, in general, the field changes 
its transverse distribution as it travels through the waveguide because dif- 
ferent modes travel with different propagation constants and different group 
velocities. The following expansion of the propagation constant k n , 




kn = \\kl - ^- = A^l - -^ - h (1 - --^j , (13) 

shows, however, that within this approximation for z = ^o^_ ^ e j n iti a l field 
is obtained. Obviously, requirements for above linear approximation are not 
met for an arbitrary k n . Higher modes have to be prevented from contributing 
to the image, because only then approximation of the square with accuracy 
to the linear term is sufficient. The reason why revivals appear is that 7$ 
factors are all of the form constant (characterizing the system) times integer. 
The question arises whether in cylindrical mirror waveguides similar analytical 
formula for the Talbot distance can be obtained. 

Let us note here that necessity of taking care of paraxial approximation is the 
main difference between optical and quantum mechanical Talbot revivals. In 
the case of infinite potential well in quantum mechanics, eigenvalues of the 
system are given by rL ^- and no excluding of modes (wave functions) that 
follows from linearization of square root is needed. 



4-2 Cylindrical mirror waveguides 



Solutions of the wave equation in cylindrical mirror waveguides corresponding 
to angular dependance e ±lv<p , where v denotes a number of Bessel function of 
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a radial solution. Propagating constants are of the form of y fc 2 — ^ for TM 



modes and y k$ — ^p- for TE modes, where j vn denotes the nth root of Bessel 
function J u (x) and j' vn the nth root of its derivative and a is a waveguide 
radius. 

In this case not only a linear approximation of square root is needed but also 
approximation for zeros of Bessel function j un or j' un . Standard asymptotical 

expansions for j m , j' vn are given by j un ~ nn + (u — |J | and j' un ~ nn + 

( z/ — | j |. They are believed to be good enough for n > v (or in more rigorous 
manner for n > 2z/). Using these formulas we can repeat procedure from Eq. 
(13) and obtain: 
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(4n + 2i/- l) 2 = 8(2n 2 + 2nz/-n+ — j +1. 

Omitting common phase factor exp (ik^z + i x 2 g^ z) we find that for a given 
wavelength A and waveguide radius a at distance zt = 16a 2 /A and its in- 
teger multiples Talbot revivals are to be obtained. Similarly, for TM modes 
approximation for j' un leads to the same Talbot distance zt = 16a 2 /A. 



4-2.1 Examples of revivals in cylindrical mirror waveguides 



To present the example of the revivals of initial intensity, I org (0), at the Talbot 
distance and its multiples we have chosen intensity function of the form: 

i sin2vrp for p e [0, 0.5), y?G[0,27r[ 

■l org (U ) — \ 

± sin 2rcp for p e (0.5,1], <p e [0, 2tt[. 

Figures 7, 8 show cross-sections for arbitrary tp of this initial light intensity 
1(0) and its Talbot revivals obtained for ak = 10 3 and ak = 10 4 , respectively. 



Analytical function from Eq. (15) and its approximation by Bessel functions 
corresponding to TM 0i modes are plotted together in Figures 7. a, 8. a. The 
numerical procedure allowing to approximate a given intensity function we 



14 



10 3 ; 

0.8 



0.07133 




0.1429 



c \ 0.2 0.4 0.6 0.8 1 d 0.2 0.4 0.6 0.8 1 e \ 0.2 0.4 0.6 0.8 1 

Fig. 7. I(p) = sin(27rp) for p G [0, 0.5) and I{p) = \ sin(27rp) for p e [0.5, 1]. The 
intensity of light is plotted: a) shows initial I(p) (an origin and an approxima- 
tion given by the first 50 terms of BF series, plotted together); b) shows the first 
Talbot revival at z = zt and the initial intensity shown together; c), d), e) show 
2nd, 5th, and 10th Talbot revival, respectively. 
ak = 10 4 , , , 
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Fig. 8. I{p) = sin(27rp) for p G [0, 0.5) and I{p) = \ sin(27rp) for p E [0.5, 1]. The 
intensity of light is plotted: a) shows initial I(p) (an origin and an approxima- 
tion given by the first 50 terms of BF series, plotted together); b) shows the first 
Talbot revival at z = zt and the initial intensity shown together; c), d), e) show 
2nd, 5th, and 10th Talbot revival, respectively. 

have used here is described in the "Methods" section at the end of the paper. 
Quality of this approximation, /q , is calculated using measure from Eq. 9. To 
indicate that the infidelity of approximation is measured, we shall write /q , 
while the infidelities of revivals defined in a similar way are denoted without 
this " a " superscript. 

Figures 7.b, 8.b present the first Talbot revival of the initial fields, Figures 
7.c, 8.c; 7.d, 8.d; 7.e, 8.e show 2nd, 5th, and 10th Talbot revivals, respectively. 
Above every plot the corresponding infidelities are depicted. It is very interest- 
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ing to compare plots and infidelities obtained for the same field distribution, 
but for different values of ak, as it is clearly seen that infidelities are much 
lower for smaller A to a ratio. This effect is connected with the quality of linear 
approximation of the square root from Eq. (14). The smaller percentage of all 
modes is used to construct the initial picture the lower infidelities of higher 
revivals are to be expected. 



4-2.2 Some comments on approximations used 

Revivals presented in previous paragraph are quite faithful which means that 
approximations used to predict the existence of revivals were justified. How- 
ever, please note, that the examples studied in the previous subsection had the 
following property: superposed modes were of the same angular dependence 
(they have corresponded to the Bessel functions of fixed number is). Studying 
more complicated combinations one finds out that superpositions of modes 
with different azimuthal mode numbers do not revive at the Talbot distance 
[21]. 

Numerical simulations of field propagation show that the situation is really 
interesting. As we have already mentioned, superposition of TM or TE modes 
having the same angular dependance revive quite faithfully at Talbot distance 
and only for the superpositions of "mismatched" modes revivals are not ob- 
tained. Explanation of this surprising fact is the following: the higher terms of 
the asymptotic formulae for the roots of Bessel functions, and its derivatives 
[22]: 
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are more important then we have assumed so far. In many papers and text- 
books only first terms (A, A') of above approximations are used and we have 
also limited ourselves to them in preliminary calculations but to obtain a 
faithful approximation at least two first terms should be taken into account. 
From formulas (15), (16) it is, however, clear that when higher terms are taken 
into account, finding the Talbot distance for arbitrary v and n fails, because 
different powers of n would be included. 
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So how we can observe any revivals at all? Let us take a closer look at approx- 
imation (16), i.e. one considering TE modes - keeping in mind that similar 
analysis can be made for TM modes as well. Taking two first terms of approx- 
imation (16) the following expression for j' un is obtained: 
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Ratios B' /A' and 2A'B' / A' for a wide range of parameters v and n are 
shown in Figure 9. It is seen that although the ratio B' /A' for n > v is close 
to zero and can be neglected, the percentage value of IB' jA' can be quite 
significant. 
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Fig. 9. Errors of "asymptotic" approximation of (j' un ) 2 : Notation is taken from from 
Eqs. 16, 17: Plots a) and b) show in two different scales percentage ratio of 2B' /A 1 
for v G {0, ..., 12}, n £ {1, ..., 25}. Unicolor dotted lines correspond to one value of 
azimuthal mode number v. Plots c) and d) present percentage ratio of B' /A' . It 
is seen that for n > v this term can be neglected. 

Obviously, this is the effect we were looking for as term 2A'B' = u 2 + | 
depends only on the azimuthal mode number v and not on the radial one n. 
That is the reason why modes having the same angular dependence revive 
at the Talbot distance zt = y, while for a superposition of modes with 
different angular dependence we do not obtain such revivals. In the first case 
factor exp (—i(4u 2 + 3)zr/4) is merely a phase common for all the terms of 
the superposition. When the initial H z is the superposition of a form H z (0) = 
Ju{j' un P) etUip + J ^{J'ixmP) el ^ l nen & 1 z — z t phase factors corresponding to J u 
and J M differ for v ^ \x and, consequently, they do not cancel. This relative 
phase destroys the Talbot revivals promised by less accurate approximation. 
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Similar analysis shows that superpositions of TM and TE modes will not re- 
vive at the same distance because of the relative phase. Only for TM U and 
TEqj modes this relative phase disappears and Talbot effect can be observed. 
However, due to possibility of using polarizers, restrictions imposed by neces- 
sity of choosing polarization of modes used is experimentally less demanding 
then that concerning y9-dependance. 



5 Summary 



We have discussed in details approximations that appear in the study of the 
Talbot effect in the cylindrical mirror waveguides as well as different approxi- 
mations used in the case of dielectric waveguides. We have shown that in many 
cases almost perfect revivals can be obtained and that even dephased propaga- 
tion can be used in practice. We have stressed that a phase space description 
(sometimes regarded as an unnecessary complicated representation) extracts 
an essence of the interference phenomena and that the conditions needed to 
be fulfilled to obtain Talbot revivals in a very natural way follow from a phase 
space description of interference. 



6 Methods 



The following numerical procedure was used for decomposing given initial light 
intensity I org into TM$i modes: For an arbitrary spherically symmetric TM field 
with an initial E z of a form 
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EM = YxMJMP) (18) 



we have evaluated the light intensity 1(0) and this 1(0) was expanded again in a 
Bessel- Fourier series corresponding to Jo- Numerical evaluation of some integrals 
was required at this point. In such a way general "basis" (ei, e2, ..., eso) was ob- 
tained, (every e$ being a sum of all possible pairs X{ ■ Xj with numerically calculated 
coefficients). Then, an arbitrary (^-independent intensity I org that we would like to 
propagate through the waveguide was decomposed in BF series 



50 
lorg(0) « ^CiJoiJoiP), (19) 

i=l 

and a set of quadratic equations c\ = e± , c<i = e2 , ■ ■ ■ ,050 = eso was solved numer- 
ically to find (xj)f£ 1 . Figures 7.a, 8. a can be treated as a test of faithfulness of the 



solutions founded in the procedure described above - the infidelities of approxima- 
tion /g are of the order of 1CP 2 . 
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